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Abstract 

In the high density low temperature limit, Quantum Chromody- 
namics (QCD) exhibits a transition to a phase characterized by color 
superconductivity and by energy gaps in the fermion spectra. Under 
specific circumstances the gap parameter has a crystalline pattern, 
breaking translational and rotational invariance. The corresponding 
phase is the the crystalline color superconductive phase (or LOFF 
phase). In the effective theory the fermions couple to the phonon 
arising from the breaking of rotation and translation invariance. We 
compute the parameters of the low energy effective lagrangian de- 
scribing the motion of the free phonon in the high density medium 
and derive the phonon dispersion law. 



1 Introduction 

QCD at large densities is expected to exhibit color superconductivity Jj]] from 
condensation of quark bilinears. It has recently been pointed out |2], |3|, ||, |5j, ||, 
[7|, H, |], H3 that a distinctive phase, called crystalline color superconducting 
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phase, may appear, in situations where different quark flavors have different 
chemical potentials, with potential differences lying inside a certain window. 
Crystalline color superconductivity is also expected to occur simply because 
of the quark mass differences [DJ. However in the following we shall limit 



the discussion to the case of two massless quarks with different chemical 
potentials. 

Besides its theoretical interest for the study of the phase structure of QCD 
theory, such a crystalline phase may result relevant for astrophysical dense 
systems @, [Q. A pairing ansatz similar to that of the QCD crystalline 



phase was already encountered in condensed matter (and referred to as LOFF 
phase from the name of the authors who studied it in superconductivity |L3|1 ). 
Other recent applications of the LOFF phase in organic superconductors, 
nuclear matter and cold trapped gases are respectively in Refs. |14] , Jl5| and 

In the crystalline phase, rotation and translation invariance are both 
spontaneously broken. The ensuing goldstone structure and the effective 
lagrangian were studied in [|J. For plane- wave behavior of the crystalline 
condensate the conclusion was that a single goldstone appears, the phonon. 
The appearance and description of a single goldstone in a situation where 
both translations and rotations are spontaneously broken was found in |4[] to 
originate from a functional relation to be satisfied by fictitious fields associ- 
ated to broken rotations such as to assure small long wavelength fluctuations 
of the order parameters. 

The main purpose of this paper is to calculate at first order in the deriva- 
tive expansion the parameters of the effective lagrangian for the crystalline 
phase and derive the phonon dispersion law. Due to the fact that quark 
condensation singles out one direction in space, the phonon dispersion law is 
anisotropic. 

The paper is organized as follows. In Section |^ we review the velocity- 
dependent effective lagrangian approach we shall use in this paper. In Section 
[| we review the basic facts concerning the crystalline super-conducting phase 
and the modifications that are needed in the effective lagrangian to deal with 
the anisotropy induced by the LOFF condensates. In Section |] we derive the 
gap equations for the BCS and the LOFF condensates by a Schwinger Dyson 
equation derived in the framework of the effective theory. In Section ^ we 
introduce the phonon field associated with the breaking of rotational and 
translational invariance and we write down its coupling to the fermion quasi- 
particles. In Section |] we derive the parameters of the effective lagrangian 
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for the phonon. Finally in Section ^ we draw our conclusions. 



2 Review of the velocity-dependent effective 
lagrangian 

We shall employ in the sequel the effective lagrangian approach based on 
velocity-dependent fields. It has been developed in [T7], [T8|, [T!| for the 
2SC and CFL models and in [|J for the LOFF phase. In this Section we give 
a brief review of it. Let us consider the fermion field 

and decompose the fermion momentum as 

p» = ^ + (2) 

with f M = (0,v) and v the Fermi velocity (for massless fermions \v\ = 1); 
finally is a residual momentum. 

By the decomposition (0) only the positive energy component ip + of the 
fermion field survives in the lagrangian in the \x — > oo limit while the the 
negative energy component ^_ can be integrated out. These effective fields 
are velocity dependent and are related to the original fields by 

1>(x) = ^e^[il> + (x)+1>-(x)], ( 3 ) 



where 



1 ± a ■ v f di r +OQ d£ _ itx 



Mx) = —j— L ^ e -** M e). (4) 

Here means an average over the Fermi velocities and 

(7 

ij)±(x) = ip±,$(x) (5) 

are velocity-dependent fields. The positive energy effective field carries color 
and flavor indices: ip + = ip+,i a where a = 1,2,3 is a color index and i is a 
flavor index . In this paper i — 1,2 (two flavors). The integration measure 
in the isotropic (CFL or two-flavor 2SC models) is as follows 

t df r dfi r+ 5 r 27T r+ x r+ 5 

J J 5 Ml1 = frCLCfI (2?r) Jo ^ I 1 dC ° S6 J 8 d£l1 (6) 



(2tt) 3 J (2tt 
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where In is along v and d£± is a surface element, which on a sphere of radius 
/i is equal to /i 2 dfl = fi 2 d(f>d cos 9; £\\ is limited by \£\\\ < 5, where 5 is such 
that (A the BCS gap) 

A < 5 < n . (7) 

The condition \£»\ < 5 is equivalent to the requirement 

\\p\-»\<5 (8) 

once one recognizes that, due to the spherical symmetry, one can always 
choose £ in the direction of v. On the other hand the condition ([?]) has a 
twofold implication: 5 < means that one is considering only degrees of 
freedom near the Fermi surface, while, by the choice A C 5, the physical 
quantities become independent of the cut-off procedure. 

In the case of anisotropic superconductivity, such as the LOFF phase, the 
situation is different and will be discussed in the next section. 

It is useful to use a different basis for the fermion fields by writing 

^P+,ia = 2^ — 7^-^+ 0> " = 1, 2) 

A=0 V2 

^+,13 = <f% 

^+,23 = , (9) 

where a a are the Pauli matrices for A = 1, 2, 3 and er = 1. (p+ are positive 
energy, velocity dependent fields: 

<4 = <* • ( 10 ) 
and represent left-handed Weyl fermions. In view of the introduction of a 
gap term (see next section) we introduce also the fields 

= • (ii) 

tp± should not be confused with the positive and negative energy fields; they 
are both positive energy fields, but are relative to opposite velocities. 

The effective lagrangian Co for gapless fermion fields at high density can 
therefore be written as follows []: 

£o = EE vtfdV ■ d) + (L - R) . (12) 

v A=0 



^ere V = (1, v) while = (1, -iJ). The derivation of is in Q. It can be 
obtained noting that, in the \x —> oo limit, terms bilinear in the fermion fields must have the 
same velocity by the Riemann-Lebesgue lemma, while the spin matrices are substituted 
by the vector by making use of simple algebraic identities |L7| . 



4 



Using the fact that the average over velocities is symmetric we write: 



1 5 

■o = o E E (v- 



•i f (iy • d) + ^(iv ■ d) <pt) . 



(13) 



Introducing now 



1 ( < 



) 



(14) 



the lagrangian can be written as follows: 



£ = EE^f* 

v A=0 \ 





iV-d 




(15) 



To this lagrangian involving only the left handed fields, one should add a 
similar one containing the right handed fermionic fields. Also notice that 
we have embodied the factor 1/2 appearing in Eq. (|13|) in the definition of 
the sum over the Fermi velocities. This extra factor cannot be eliminated 
by a redefinition of the fields if>± since it corresponds to a genuine doubling 
of the degrees of freedom; its appearance is needed because, by introducing 
the fields with opposite v, if)- = if>+,-v we must integrate only over half solid 
anglef]. 

3 Crystalline colour superconductive phase 

As shown in when the two fermions have different chemical potentials fi\ ^ 
/12, for Sfj, of the order of the gapQ the vacuum state is characterized by a non 
vanishing expectation value of a quark bilinear which breaks translational and 
rotational invariance. The appearance of this condensate is a consequence 
of the fact that for ^ and in a given range of 5fj, = — //2I, the 
formation of a Cooper pair with a total momentum 



2 In [[19) another factor 1/2 appears in the lagrangian Cd\ here we simplify the notation 
and get rid of it by a redefinition of the fermion fields. 

3 The actual value of the range compatible with the presence of the LOFF state depends 
on the calculation by which the crystalline colour state is computed. Assuming a local 
interaction as inB produces a rather small interval, which is enlarged by assuming gluon 
exchange, as in M. 



Pi + p 2 = 2q ^ 



(16) 
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is energetically favored in comparison with the normal BCS state. The anal- 
ysis of shows that two different, but related condensates are possible, one 
with the two quarks in a spin zero state (scalar condensate) and the other 
one characterized by total spin 1 (vector condensate). In the BCS state the 
quarks forming the Cooper pair have necessarily S = 0; as a matter of fact, 
since the quarks have opposite momenta and equal helicities, they must be 
in an antisymmetric spin state. This is not true if the total momentum is 
not zero and the two quarks can have both S = and S = 1. 

The possible form of these condensates is discussed in p|; these authors 
assume only two flavors and make the ansatz of a plane wave behavior for 
both condensates: oc e 2iq ' x . Though more complicated structures are possible 
we will make the same hypotheses. The results of B £1X6 cLS follows: 

- < 0\e tl e af33 ^ a (x)C^(x)\0 >= 2T L A e 2 ^ ; (17) 



besides the scalar condensate ([17]) one may have a spin 1 vector condensate, 
symmetric in flavor: 

i < 0\4e a ^ a (x)Ca 03 ^(x)\0 >= 2Y L B e 2 ^ . (18) 

The effect of the two non vanishing vacuum expectation values can be taken 
into account by adding to the lagrangian the term: 

r — r( s ) _i_ — 

J~>A — *-A ~r ^"A 



2 . 

._ e «/33 



1&(x)C (AW + a ■ nA^4) ^(x) + h.c. , (19) 

which includes both the scalar and the vector condensate. 

Let us consider the lagrangian term relative to the scalar condensate. We 
shall write it as follows: 

4? = e 2l ^e^\ l3 ^ T ia {x)C^{x) -(L^R) + h.c. . (20) 

Here ip{x) are positive energy left-handed fermion fields discussed in the 
previous section. We neglect the negative energy states, consistently with 
the assumption discussed in the previous section. 

In order to introduce velocity dependent positive energy fields 
with flavor i we have to decompose the fermion momenta according to (fj): 

Pj = m + ?j (i = i,2). (2i) 
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Therefore we have : 

= 2~ X! exp{ix- a(vi, v 2 , q^}e ij e a ^ii +> -^.i a {x)C^ +t -g j .jp(x) 

V1,V2 

-{L -> i?) + h.c. , (22) 



where 



a(#i, u 2 , g) = 2g - /ii^i - /j 2 v 2 ■ (23) 



We choose the z— axis along g and the vectors v\, v 2 in the x — z plane; if 
«i, «2 are the angles formed by the vectors v±, v 2 with the z— axis we have 

oc x = — /ii sin«i cos0i — /i2 sin a 2 cos02 , 
Oi x = — /ii sinai sin0i — /i 2 sin«2 sin0 2 , 
ot z = 2q — /ixcosai — /i 2 cosa 2 . (24) 

Because of the already quoted Riemann-Lebesgue lemma, in the //i, \i 2 — > oo 
limit the only non vanishing terms in the sum correspond to the condition 



ot x = ot y = . (25) 



Let us introduce 



2 



8v = -^=^. (26) 



ai = a 2 + 7T + - , (27) 




Eq. fl25|) implies 



i.e. the two velocities are opposite in this limit. 

The condition on a z depends on the behavior of q in the limit \x — ■> oo. 
If we might take the limit q ~ — > oo as well, we would also have ot z ~ 0; 
we will make this approximation which is justified by previous analysis 0, 
which points to a value of q ~ 1.25/x and 5/x ~ 0.7 Abcs ^> A^' - "). These 
conditions have to be justified also in our approach and we will do it in the 
next section. Notice that this approximation considerably simplifies the loop 
calculation in presence of the LOFF condensate, which is usually performed 
by defining a blocking region, where the LOFF condensate is impossible, and 
an allowed, pairing region, whose definition is mathematically involved. This 
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approximation may be therefore useful in complex calculations such as the 
present one, even though it can produce numerical results somehow different 
by the more precise calculations given in ||. We have: 

5fi q 2 - Sa 2 
cosai = h 



q qn 

Sa q 2 — 5a 2 , , 

cosa 2 = + — + - — , 28 

q qn 



which, together with (f25|) , has the solution 

. Sa e 

a 2 = Oq = arccos — — - , (29) 

q 2 

«i = a 2 + 7r — e (30) 

with 

e = 2 Vq2 - 6 ^ . (31) 

Therefore, as anticipated, e = O [ — ] and in the limit a — > oo the two 

W 

velocities are almost antiparallel: 

vi ~ -v 2 . (32) 

Putting 

ll>+,±Sf,ia(x) = 1p±^- ia (x) , (33) 

eq. Q2"2j) becomes: 



£a = E ^ a ^ + v; ta (x)C^ J „(x) ~(L^R) + h.c. (34) 

v 

In a similar way the term corresponding to the vector condensate in the 
lagrangian can be written as follows: 

4? = -^E^^;«WC(?-n)^ ; ^(x) -(L^R) + h.c. (35) 

v 

where n = q/\q\ is the direction corresponding to the total momentum carried 
by the Cooper pair and we have used ip^Ca ■ nip + = v ■ "nip^Cip^ . 
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The procedure of the limit q — > oo may produce a finite renormalization 
of the gap parameters. To show this let us consider in more detail how 
the second of the two equations ( |28|) can be obtained. In order to take the 
limit in the lagrangian for fi —>■ oo or q — > oo we perform a smearing in the 
space time integrations; for example, considering the z axis we introduce a 
smearing \Az\ and we write: 

e 04h* _» 1 f Z+lAZl/2 ^ qh y d _ ^Hz Bin(gfe| Ag[) _ 

|Az| Jz-\Az\/2 qh\Az\ 
= ^ hz 6 R (h) (36) 



1? 

We have put 



i? = q\Az\ (37) 



and introduced the "fat delta" 5r(x) defined by 

sin(Rx) 

Sr[x) = , (38) 

7TX 

which, for large R, gives 

5 R (x) -> . (39) 

Let us now justify the approximation (|28|). We have (02 = azimuthal angle 
of the velocity vector v 2 ) 

J d cos 9 2 e t2qh{cos 02)2 -> ^^Jdcos9 2 e i2qh(cose2)z 5 R (h(cose 2 )) 

~ 1 a iT777 -r / dcos02<5(cos0 2 - cosa 2 )(40) 

q\Az\\h'{cosa2)\ J 

In the previous equations we have defined^ 

/i(cos0 2 ) = 1 + -^ |-/i2 cos 9 2 + cos 2 J/4- -^tz-| • ( 41 ) 

In the sequel we will also use 

h± (cos 6 2 ) = 1 + ^- (^-^2 cos 2 ± sj fi 2 cos 2 2 - ASfifi^j (42) 
4 The reader should note that, by our convention, dfi < 0. 
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where one has ±1 = +sign(cos6 , 2 ). Due to eq. P0|), the condensates are 
modified as 

A - ^5 R {h ± {cos9)) A . (43) 

We handle the fat delta according to the Fermi trick in the Golden Rule; in 
expressions involving the gap parameters we make one substitution 



5 R (h±(cos6)) -> 6(h±(cos6)) 
in the numerator; correspondingly we will get a factor 



qR 



(44) 



(45) 



in the numerator (from the region cos9 < 0); the other fat delta are computed 
as follows 

7rd R (h_(cosO)) tt5 r (0) 

R > (46) 

In the basis introduced in the previous section, the effective lagrangian is 



C 



T >0 A 

5 



*-A "~ *"A 



t? A=0 



At 



i Sab V ■ d 

A A b i Sab V ■ d 



X 



B 



The matrix Aab is as follows: 

A AB = 



while for A,B = 0, 3 we have: 



(AorB = 4 or 5) 



(47) 



(48) 



7T 

Aab = -^S R (h ± (cos9) 



f A 








-Ax 


\ 





-A 


-zAi 










-A 







V Ax 








-A 


/ 



(49) 



with 



A 
Ai 



i7-nA M 



(50) 
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Let us now discuss the precise meaning of the average over velocities. As 
discussed above we embodied a factor of 1/2 in the average over velocities; 
however in the LOFF case the sum over the velocities is no longer symmetric 
and, whenever the integrand contains the matrix Aab now reads: 



2 



2ty 



if 



v = (sin 9 q cos (ft, sin 9 q sin (ft, cos 9, 



(51) 



(52) 



and 9 q given in (|29|). As discussed above, this limit corresponds to i? large. 
For future reference we observe that, in the same limit, 



v i v i = ha [ Sm /I — (SuSji + 5i 2 6 j2 ) + C ° S ® q 5i 3 5 j3 



\ 4 v J J ' 2 

V \ / 

The effective action for the fermi fields in momentum space reads: 

dH dH' 



(53) 



S = E E 

v A,B=0 



(2vr) 4 (2tt) 



-X A \nD- AB {t' A)xb{£) , 



(54) 



where £>^(f,£) is the inverse propagator, given by: 



V-£5 



.t 



AB _ ^AB 

A AB V 



5\e - £) . 



(55) 



From these equations one can derive the quark propagator, defined by 
d^' 



E 



B 



(2tt)< 



Dj B (£,£')D B c(£',£") = S AC S\£ - £") . 



(56) 



It is given by 



D AB (£,£") = (27r)W(£-£")xY^ 

c 



/ V-£5ac 

D C b(£) D cb (£) 

A A c V-£5 A c 

V D CB {£) D CB (£) J 



(57) 
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where 



d C b(() =(y-ev-e- aa*) cb 

D CB (£) = {y-iV-i- AtA) cB . (58) 

The propagator for the fields x 4 ' 5 does not contain gap mass terms and 
is given by 

D(£, t) = (2tt) 4 5\£ - t) ( {V ^ _1 ( - ) • (59) 

For the other fields x , A = 0, • • • , 3, it is useful to go to a representation 
where AA' and A^A are diagonal. It is accomplished by performing a unitary 
transformation which transforms the basis \ A into the new basis \ A defined 
by 

X A = RabX B , (60) 

with 



Rab 



In the new basis we have 





( 1 





1 \ 


1 





1 -% 










+i -1 







V 1 





-1 / 



(AAt)^ = X A 5 AB 



(61) 



(A f A) AB = ^aSab (62) 

where 

Aa = (^5 R (h ± (cose) ) 2 ((A + A x ) 2 , (A - A^ 2 , (A + A x ) 2 , (A - A : ) 2 ) 

Aa = {^ R {h ± {cos9) ) 2 ((A - AO 2 , (A - A x ) 2 , (A + A^ 2 , (A + A : ) 2 ) 

(63) 

For further reference we also define 

He = (A + A ls A 1 - A , A + A 1; A x - A ) . (64) 
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Let us finally discuss the integration limits in the residual momentum £ 
in the LOFF case. We write: 

e = £ x + £\\ , (65) 

where in is parallel to v% ~ — v% and £± is a 2-dimensional vector orthogonal 
to v 2 . Clearly V ■ £ = £ — in and V ■ I — £q + in, therefore the propagator 
only depends on £ and £\\; however, since the integration over velocities is 
anisotropic, the theory, differently from the CFL and 2SC cases, is not 2- 
dimensional, but (2+l)-dimensional. One has for the integration measure, 
whenever the integrand is proportional to the gap 

r dH r d£ ± r+* /•+«> d£o 

J (2?r) 4 ~ J (2tt) 3 J -s 11 J-oc 2^ ~ 
47T// 2 k R r+ s r+°°d£ 

= j2^YLs di "L 2^' (66) 

which implements the condition (p9|). 

Before closing this Section we note some differences between the present 
calculation and the one presented in Q. In 0] we did not work in the 
simplifying approximation a z = (see eq. (p^), which implies that the 
effective lagrangian and the fermion propagator in ||J] contain both 2q and 
S n and are therefore more involved (there is also a numerical mistake in 0] 
that renders hermitean the gap matrix Aab, which is not true. The correct 
expression is eq. (|49|)). 



4 Gap equation 

To derive the gap equation in the CFL case we write a truncated Schwinger- 
Dyson equation, similarly to the approach followed by [H. We assume a 
fictitious gluon propagator given by 

^=2^, (67) 

which corresponds to a local four-fermion coupling. Using the Feynman rules 
that can be derived by the effective theory one obtains, in the limit of eq. (p9|) 
two equations: 
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/I 2 

Al = ~ l 24A%3 ^ o/( ^ /io - > » ( 68 ) 

where 

/(A) = / = -i2vrarcsinh ( 4r I • (69) 

We fix A by using the same equation and the same fictitious gluon prop- 
agator at n — 0, where the effect of the order parameter, in this case the 
chiral condensate, is to produce an effective mass M for the light quark. This 
assumption gives the equation 

rK r? 



V*l dv JWTW- (70) 



3A% 2 Jo " VP 2 + M 2 ' 

For M = 400 MeV and the cutoff K = 800 MeV (700 MeV) we get A = 181 
MeV (154 MeV). 

On the other hand we can consider the analogous equation for the 2SC 
model which is given by 

' '(A) = ^arcsinh(^| (71) 



3A 2 7T 3 v ' 3A 2 7T 

that can be solved explicitly: 



sinh 



2/x 2 J 

With the same values of A (A = 181 MeV and 154 MeV) we find the val- 
ues for the gap parameter of the 2SC model reported in the following table, 
where we have taken 8 of the order of /i. 



5 = K- /i 


A 


400 MeV 
300 MeV 


39 MeV 
68 MeV 



Table 1. Value of the gap parameter A in the 2SC model for two values of 
the cutoff K . The chemical potential has the value jj, = 400 MeV. 
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Similar results have been obtained by with a procedure which differs 
from the present one for two reasons: They include the 0(1/ fi) corrections 
to the gap equation and use a smooth cutoff instead of the sharp cutoff K 
used here. 

The coupled equations for the LOFF case fl68f) have a non trivial solution: 



A o = A7T2-2V > A i = • ( 73 ) 



sinh 



3AV 
H 2 kn 



This is the result for the gap in the LOFF phase in the approximation 
fx — > oo, q — > oo . It is known that the crystalline phase holds in a well 
defined window of values \5fi\: \5fi\ G (Sfii, <5/i 2 ). The determination of 
is obtained by comparing the free energies of the normal and 2SC states and 
is completely analogous to the one in [Q], based on eq. (4.3) of this paper; 
one gets, in the weak coupling limit 8fi\ ~ 0.71A25C- On the other hand 
S/j,2 is the limiting value of \5fi\ between the LOFF and the normal phase; 
at \5[i\ ~ 5fi2 one expects a second order phase transition and Ao — > 0. 
In our approximation the phase transition occurs because, as it is evident 
from (|28|), for sufficiently large 5/i and fixed q (in the analysis of 0, in the 



LOFF window q is constant m 0.90A 2 sc* ~ 35.3 MeV) cos 6^ is not inside 
the integration region in cos# and the only solution of the gap equation is 
A = 0. Therefore, within this approximation, <5/i 2 ~ q ~ 0.90A 2 5C) which 
compares favorably with the result ~ 0.754A 2 5c of 0. 

Going beyond this approximation implies the use of the eq. ( p5l) alone, 



instead of ( [25] ) and (P5]); the result ( PB] ) is indeed more robust, based as it 
is only on the limit /i — > 00 and not also on q — > 00. Now the gap equation 
would assume the form (5 = fi): 



1 



( J dz 6ji(h_(z)) arcsinh 



3A 2 7T 2 i? 

+1 



l6 R (h-(cos9) A 



/ dz8 R (h+(z)) arcsinh ——7 — — — ) • (74) 

xd R (h + (cos6) A 







Also for this approximation we find non vanishing gaps; for example at R 1 
and for SfM = 5fii, we get A ~ 5.5 MeV; in this case 8^2 ~ l-25A 2 5c! for 
R ~ 1.5 we find smaller values for the gap: Ao ~ 0.25 MeV and 5/i 2 ~ 
0.83A 2 5<7. We observe explicitly that in the limit of large R one gets, as a 
solution of (|74|) , the eq. (|73|) with 5 = \i. 
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We finally notice that we find Ai = 0, whereas in the more precise calcula- 
tion of ref. it was found a value A x 7^ ; but very suppressed with respect 
to Ao- Presumably the difference comes from our expansion at the leading 
order in /1, and this would explain the suppression of the vector condensate. 

Notice that, even though our analysis leads to Ai = 0, we will consider 
in the sequel for completeness the general case where both condensates are 
present. 



5 Phonon-quark interaction 

Let us consider again the breaking terms in (|T9"D, i.e. and Cj±- This 
lagrangian explicitly breaks rotations and translations as it induces a lattice 
structure given by parallel planes perpendicular to n: 

n-x = 7 ^- (k = 0, ±1, ±2, ...) . (75) 

We can give the following physical picture of the lattice structure of the LOFF 
phase: Due to the interaction with the medium, the Majorana masses of the 
red and green quarks oscillate in the direction n, reaching on subsequent 
planes maxima and minima. The lattice planes can fluctuate as follows: 

n ■ x — > n ■ x + -— , (76) 
*Qj 

which defines the phonon field, i.e. the Nambu-Goldstone boson that is 
associated to the breaking of the space symmetries and has zero vacuum 
expectation value: 

(0)o = O. (77) 



Let us observe that, because of (|75f ), we have a fluctuating field 4>k for any 
nodal plane, where 

fc = <f)(t, x, y, z k = kn/ q) (78) 
Also the vector n can fluctuate: 

n -> R, (79) 

where the field vector R satisfies 

\R\ = 1 , (80) 
(R)o = ft. (81) 
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Also the vector fields R represent, similarly to (|75|), a collection of fields: 



R = R k . (82) 

However we are interested to give an an effective description of the fields 
<pk and Rk in the low energy limit, i.e. for wavelengths much longer than 
the lattice spacing ~ 1/g; in this limit the fields <pk and Rk vary almost 
continuously and can be imagined as continuous functions of three space 
variables x, y and z. Therefore we shall use in the sequel the continuous 
notation and R and postpone a discussion on this aspect to the end of the 
present Section and to the subsequent Section. 

The vector field R behaves as a vector under rotations while being invari- 
ant under translations. If we take the z-axis pointing along the direction n, 
we can write 

(R)i = iM£i,6) = (e^ Ll +^)/^) j3 (83) 
with L the generators of the rotation group in the spin 1 representation, i.e. 

(Li).,-* = -ie ijk (84) 

and £i, £2 are fields. At the lowest order in these fields we find 

R^n + 5n (85) 

with 

Sn= +f ,0) . (86) 

V JR JR J 

Since the fields Ri transforms as a vector under rotations, it follows that the 
rotational symmetry is restored by the substitution (|79|). 

If the direction q appearing in the exponent in (|H]) were different from 
the direction n of the spin 1 condensate in the same equation, the physical 
situation would not differ from magnetic materials; we would need a NGB 
to take into account the fluctuations of the scalar condensate (|I7]), i.e. 

A e 2i ^ (87) 

and other independent fields 5n, behaving as spin waves, i.e. making a pre- 
cession motion around the direction n. Here, however, the situation is differ- 
ent, because the two directions coincide. Therefore the value of R is strictly 
related to the symmetry breaking induced by the scalar condensate. 
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We have shown in 0] that in order to describe the spontaneous breaking of 



space time symmetries induced by the condensate (|8~7 ) one NGB is sufficient. 
As we stressed in [|J , rotations and translations at least locally cannot be seen 
as transformations breaking the symmetries of the theory in an independent 
way, because the result of a translation plus a rotation, at least locally, can 
be made equivalent to a pure translation (]. 

As a matter of fact at the lowest order we write 

^ T 

{n + Sn) ■ (x + 5x) — n ■ x R ■ x + n ■ x , (88) 



2qf " ' v ' 2q 

where we have introduced, as in Q, the auxiliary function T, given, in the 
present approximation, by 

T = 2qn-5x . (89) 

Now the lattice fluctuations must be small, as the phonon is a long- 

2qf 



wavelength small amplitude fluctuation of the order parameter. From ([88]) it 
follows that T must depend functionally on R, i.e. T = which, using 

again (§3), means that 

^n.t+±- r R-x + ^ = G[R}. (90) 

The solution of this functional relation has the form 

R = h[<$>] (91) 

— * 

where h is a vector built out of the scalar function $. By this function one 
can only [] form the vector V$; therefore we get 

^ = ^' (92) 



5 Recently the problem of the Goldstone theorem for Lorentz invariant theories with 
spontaneous breaking of space-time symmetries has been considered in [ p2| . 

6 In principle R could depend linearly on a second vector, x, but this possibility is 
excluded for the assumed transformation properties under translations. 
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which satisfies fl80|) and, owing to (|77|), also (PT|). In terms of the phonon 
field (f) the vector field R is given (up to the second order terms in </>) by the 
expression 



R = n + 



2fq 



V0 — n(n ■ V0) 



+ 



71 



8f 2 q 2 



3(n-V4>) 2 - \V(f>\' 



Af 2 q 2 



[n-V<p) . 
(93) 

Note that at the lowest order the fields £i, £2 introduced in fl36|) are given by 



+ 



2fq dx 
1 d0 

2fqdy 



(94) 



We observe that these relations do not involve d<fi/dz. In fact we are defining 

7T 

the phonons as fluctuations of the lattice planes z = k—. In the case of a 

q 

single plane it is clear that the phonon propagates only along the plane (say, 
x and y directions), and the only motion along z would be a translation 
corresponding to a zero mode. However in our problem we have infinite 
Goldstone bosons propagating on the planes z = z\. = kn/q. That is we have 
a collection of fields as in eq. (|78|). The presence of many planes gives rise to 
a finite energy mode corresponding to the variation of the relative distance 
among the planes. Phrased in a different way, different phonon fields interact 
each other since the fermions propagate in the whole three-dimensional space. 
Translating the fields in momentum space we have to take a Fourier series 
with respect to the discrete coordinate z^. This introduces in the theory a 
quasi- momentum p z through the combination exp(ip z Zk)- Therefore p z and 
p z + 2q define the same physical momentum and we can restrict p z to the first 
Brillouin zone — q < p z < q. As already stressed, this point will not affect 
our final discussion since we will be interested to energies and momenta much 
smaller than the gap. 

The interaction term with the NGB field is contained in 



C 



V 

(95) 

Notice that we have neglected the breaking of the color symmetry, which has 



been considered elsewhere p3 |. At the first order in the fields one gets the 
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following three-linear coupling: 

J y 

1 



Y^v- V<P-n{n-V<P) A^a^^C^ 



2/9 

-(L^R) + h.c. 
We also write down the quadrilinear coupling: 



(96) 



C 



J y 



J 1, 

1 



E 



n(n • V0)J A^a^^C^^ 
2 ( 3 (n-V0) 2 -|V0| 2 )-(^-V0)(n-V0) 



v ■ n 



4/V ^ 

xA^a^^C^^- (L ^ R) + h.c. 



x 



(97) 



We notice that in (p5|), (|96"l) and (|97|) the average over the velocities is defined 
according to (|oT|). 

It is now straightforward to re-write these couplings in the basis of the x 
fields. One gets 



a + u = E E x 

v A=0 



93 + 9a 



X 



(98) 



Here 



93 
94 



i(pA 



AB 



f 



?A 



AB 



2P 



<tabO[4>) 



+ CT AB \jOW\ + QWi 



(99) 



with 



6[<j>] 



2fq 
AM 

4/V 



V0 - n(n • V0)] A (v) 

'3(tT. V0) 2 - |V0| 2 j -(i?- V0)(n- V0) 
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,(100) 



(Tab 



( -1 \ 

-i 

+i 

V +1 / 



(101) 



For the following it is important to observe that the term in g 3 and g± pro- 
portional to Aab arise from the expansion of expi(p/f alone, whereas the 
terms proportional to <jab get also contribution from the expansion of R in 
the vector condensate. 

Before concluding this section, let us come back to the lattice structure 
given by the parallel planes of Eq. (|75|). The effective action for the field 4>, 
S[(j)], is obtained by the lagrangian as follows 



, 7T + °° 

S= / dtdxdy — £(<f)(t,x,y,kiT/q) , 



(102) 



In the action bilinear terms of the type c^k'Pk' with k ^ kl may arise. In the 
continuum limit this terms would correspond to derivatives with respect to 
the z direction 

We can use the same lagrangian £3 + £4 given in (|98|) with the integration 
J d A x provided we multiply the r.h.s. of (^) by the factor 



measure 



t(z) 




(103) 



As already observed the phonon momentum is a quasi-momentum, i.e. p z 
and p z + 2q correspond to the same z— component of the momentum. Due 
to this remark one can make use in the calculations of the formula 



E 



Jkirp/q 



(104) 



6 Effective lagrangian for the phonon field 

To introduce formally the NGB in the theory one can use the gradient expan- 
sion (see e.g. p4j), i.e. a bosonization procedure similar to the one employed 
in to describe the goldstones associated to the breaking of 5*^7(3) in the 
CFL model. First one introduces an external field with the same quantum 
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a) b) 
Fig. 1. Self-energy (a) and tadpole (b) diagrams. 

numbers of the NGB and then performs a derivative expansion of the gener- 
ating functional. This gives rise to the effective action for the NGB. At the 
lowest order one has to consider the diagrams in Fig. 1, i.e. the self-energy, 
Fig. la, and the tadpole, Fig. lb, whose result we name U(p) s . e . and H(p)tad 
respectively. 

To perform the calculation one employs the propagator given in Eq. (|57|) 
and the interaction vertices in (|98|) . The result of the calculation of the two 
diagrams at the second order in the momentum expansion is as follows: 

^ 2 fj2 B \ 4A c 4\ c v-ev-e 



4A c nam " hr w ® 



where 
\i c defined in (BJ1) and 



Dm V Q J \DW) D c (t), 

x - P 2 y + 2pI ~ 2p • tfp*) ] (105) 

A?(*) = *(! - ^ - A c + «e , (106) 



w(p) = p ■ v — (p ■ n)(v ■ n) . (107) 



One can easily control that the Goldstone theorem is satisfied and the phonon 
is massless because one has 

n(o) = n(o) s . e . + n(o) iad = o . (ios) 
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In performing the integration we take into account the circumstance that 
^ ^ At the second order in the momentum expansion one gets 



n(p) = - 



27rV : 



Here 



UJ' 



\p) ( 2 ^2 arcsinh—^— 

V 2 c=0 |/i c | 



+ 



+ 



AM 3 5 

— X! x arcsmri — 



2g 2 



2 



AM 

1 — log J (c<j 2 (p) + {7 • n$(p) 



where we have taken A^ <C A and 

^(p) = (3Pz — P 2 ) v ■ n — 2p- vp z 
After averaging over the fermi velocities we obtain 



n(p) = - 



2 2 / 2 i 2\ 2 ' 

P - V ±(Px+Py) ~ V \\P. 



2 2 

"1 



4tt 2 / 2 

One obtains canonical normalization for the kinetic term provided 



H 2 k 



R 



2tt 2 



On the other hand 



l - sin 2 e q + (l - 3 cos 2 9 q ) (l - log |^ ( 



aw 



= cos 6 n 



(109) 



110) 



111^ 



(112) 



(113) 



(114) 
(115) 



We note that f 2 depends on k,R. We expect that kn is in the range (1,2), 
where 1 corresponds to a smearing distance equal to the lattice spacing 
|Az| = ir/q, and 2 corresponds to realistic values for the gap. It is re- 
markable however that the phonon dispersion relation is not affected by this 
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uncertainty. We also notice that the only source of anisotropy is in the 
asymmetric average over the Fermi velocity, that is the integration over a 
circle instead of the integration over a sphere. In fact, if we were going 
to average the previous result over cos 6 q we would get the isotropic result 
v j_ — v 2 — 1/3. In particular the contribution from A^ 2 averages to zero. 
In fact, this contribution is just the one arising from the terms in g 3 and 
proportional to ctab- For the sake of the argument let us forget for a mo- 
ment the relation between $ and R. In this case the previous contribution 
is nothing but the "zero-momentum" contribution to the two point function 
for the fields £j (or Sn) of eq. fl83|). In the actual calculation this is a second 
order term in the momenta simply because the fields ^ are derivatives of 4>. 
It is a simple consequence of the symmetries to see that, after averaging the 
velocity over all the Fermi sphere, this term is necessarily proportional to 
£ 2 + £ 2 pa (5n) 2 . In fact, the only vectors involved are 5n, n and v. Since 
5n ■ n = 0, after averaging the result is necessarily proportional to (5n) 2 . 
However, the Goldstone theorem tells us that Goldstone fields should have 
no mass. This is true also for 8n, since being a gradient of it still satisfies 
the same wave equation of (p. This explains the particular expression for this 
term, which has to average to zero integrating the velocity over all the Fermi 
sphere. 

In any event the term proportional to the vector condensate A^ is very 
small and the result for v\ can be approximated as follows 

wi~^sin 2 ^. (116) 
In conclusion, the dispersion law for the phonon is 

E(p) = ^vl(pl+pl)+vf lP l (117) 

which is anisotropic. Besides the anisotropy related to v± ^ v\\, there is 
another source of anisotropy, due to the fact that p z , the component of the 
momentum perpendicular to the planes ([75]), differently from p x and p y is a 
quasi momentum and not a real momentum. The difference can be better 
appreciated in coordinate space, where the effective lagrangian reads 




However, in the long distance limit £ » ir/q, the set of fields <pk(x,y) 
becomes a function <p(x, y, z) and the last term can be approximated by 
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7 Conclusions 



The formalism we have used employs velocity dependent fields and it has 
been adapted to describe anisotropic phases, such as the crystalline phase of 
interest here. We have considered QCD with two massless flavors. For the 
two-quark condensates ( scalar and vector) we have made the currently used 
plane-wave ansatz. For different values of the chemical potentials of up- and 
down-quark the two fermion momenta require independent decompositions. 
The quark propagators can be derived from the effective action for the fermi 
fields, as provided by the velocity dependent formalism, neglecting antipar- 
ticles (with respect to the Fermi sphere). A careful discussion is required 
for the integration over velocities, due to the anisotropy of the LOFF phase, 
implying that the effective theory is (2+l)-dimensional, differently from the 
situation that applied to CFL or to 2SC where it was (l+l)-dimensional. 

We have derived the gap equation for the normal BCS state and the 
LOFF phase within the same formalism. Our approach is based on a few 
approximations, most notably the neglecting of the negative energy states 
and sub leading terms in the /J, — > oo limit. These approximations consid- 
erably simplify the formalism and render it suitable for more complicated 
computations such as those arising from crystalline patterns more involved 
than the simple plane wave ansatz assumed in this paper. In spite of this 
approximations we are able to obtain results that qualitatively agree with 
more complete calculations existing in the literature. 

The breaking of rotations and translations induces within the plane-wave 
ansatz a lattice structure of parallel planes corresponding to maxima and 
minima of the fermion Majorana masses. Fluctuations of the lattice planes 
define the phonon field. The phonon momentum has to be regarded as a 
quasi-momentum, but one can restrict to the first Brillouin zone. The formal 
introduction of the Nambu-Goldstone boson has been performed through 
bosonization, leading to the effective boson action. 

The main result of the present paper is the application of the formalism to 
derive the parameters of the phonon field, i.e. the goldstone field associated 
to the breaking of the rotational and translational invariance in the crystalline 
color superconductive phase of QCD. Lowest order calculation requires self- 
energy and tadpole diagrams. They have to be calculated in terms of the 
fermion propagators and the three-linear and quadrilinear phonon-quark in- 
teractions. Besides explicit verification of the zero value of the phonon mass, 
we have derived the values of the phonon constant / and the expressions 
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for the transverse and parallel velocities. The resulting phonon dispersion 
relation is anisotropic, both because of the difference in the velocity compo- 
nents, and of the quasi-momentum aspect of one momentum component. In 
coordinate space, the effective lagrangian exhibits explicitly such anisotropy 
in a most apparent form. 
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